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Lynds 1622: a nearby star forming cloud projected on 
Orion B? 
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ABSTRACT 

We present results of optical spectroscopic and photometric observation of the prc- 
main sequence stars associated with the cometary shaped dark cloud Lynds 1622, and 
12 CO and 13 CO observations of the cloud. We determined the effective temperatures 
and luminosities of 14 pre-main sequence stars associated with the cloud from their 
positions in the Hertzsprung-Russell diagram, as well as constructed their spectral 
energy distributions using optical, 2MASS and Spitzer IRAC and MIPS data. We 
derived physical parameters of L 1622 from the molecular observations. Our results 
are not compatible with the assumption that L 1622 lies on the near side of the Orion- 
Eridanus loop, but suggest that L 1622 is as distant as Orion B. At a distance of 400 pc 
the mass of the cloud, derived from our 12 CO data, is 1100 Mq, its star formation 
efficiency is ~ 1.8%, and the average age of its low-mass pre-main sequence star 
population is about 1 million years. 

Key words: ISM: clouds; ISM: individual objects: L 1622; stars: formation; stars: 
pre-main sequence 
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1 INTRODUCTION 

Studies of the large-scale properties of star forming regions 
are important for understanding the interstellar processes 
governing star formation. Census of young stellar objects 
(YSOs) born in a cloud, their mass, age and space distribu- 
tion, together with the mass of the star forming cloud are 
key data for assessment of the efficiency and time scale of 
star formation. 

The Orion OBI association and its associated molecu- 
lar clouds represent a distinguished target for star forma- 
tion studies, being the nearest giant star forming region 
where various stages of the star forming processes across 
the whole stellar mass spectrum, as well as interactions be- 
tween the consecutive generations of stars can be studied in 
detail. The target of the present paper, Lynds 1622 and the 
neighbouring Lynds 1621 are small dark clouds at the low 
Galactic latitude edge of the giant molecular cloud Orion B. 
These two clouds together are termed as Orion East by 
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iHerbig fc Raol(|l972l h who discovered five Ha emission stars 
associated with L 1622. Maddalena ct al. ( 1986) found that 
the radial velocity of this cloud is about ulsr = +1 kms -1 , 
contrary to the characteristic radial velocity of +10 kms -1 
of Orion B. The asymmetric, cometary shape and the bright , 
ionized rim apparent in the SHASSA (jGaustad et al.ll200ll ) 
image of L 1622 suggest its inter action with the hot stars of 
Ori OBl.[C asassus et al.l (|2006t ) present centimeter contin- 
uum image of L 1622, and attribute the radiation to spinning 
dust grains in a limb-brightened shell of the cloud, where the 
incident UV radiation from the Ori OBlb association heats 
and charges the grains. 

Low-mass star formation in L 1622 is indicated by sev- 
eral associated T Tauri stars and candidates. The reflec- 
tion nebula VDB 62, illuminated by the K2 type w eak-line 
T Tauri star HBC 515 (HP 288313, V1793 Ori) (jRacind 
ll968l : lHerbig fc Belll[l98sh can be found near the bright rim 
of the cloud. Spectroscopic and p hotometric variability o f 
this star were studied in detail by iMekkaden et al.l d2007l), 
and it was found to be a triple system by Reipurth et al.l 
(2008a). In addition to this WTTS, there are six known clas- 
sical T Tauri stars (CTTS) in L 1622: HBC 188 (LkHa 334), 



HBC 189 (LkHa 335), HBC 190 (LkHa 336), HBC 516 
(Lk Ha 336c), HBC 191 ( LkHa 337), and HBC 517 (L 1622- 
15) jHerbig fc BelJll988h . The L kHa 336 triple system was 
studied bv lCorreia et all (|2006l ). IOgura fc Hasegawal (ll983l ) 
detected 16 Ha emission stars (L 1622-1 to L 1622-16) in an 
objective prism survey over the surface of L 1622, including 
the above HBC stars except LkHa 336c and HD 288313. The 
pre-main sequence (PMS) nature of 11 objects has not yet 
been confirmed. 

A further signpost of low- mass star formation is th e 
Herbig-Haro object HH 122 (|Reipurth fc Madsenl [l989h . 
wh ose probable exciting so u rce, HH 122 VLA 1 w a s dete cted 
by iRodriguez fc Reipurthl (| 19941 ). iLee fc Mversl (|l999h re- 
port four optically selected cores of the L 1621 / L 1622 
region: L 1621-1, L 1621-2, L 1622 A and L 1622 B. 
L 1622 B is associated wit h the optically invi si ble s ource 
IRAS 05522+0146, while ILee. Myers. Tafallal (l200ll) list 
L 1622 A as an infall candidate. Recently iReipurth et alJ 
(|2008bl l identified 28 young stellar objects (YSOs) in the 
Spitzer IRAC images of L 1622, 10 of which they classified 
as protostars. 

The velocity difference between Orion B and L 1622 
raises the question whether both clouds are located in the 
same volume of space, or L 1622 is a foreground object, un- 
rela ted to the Orion B mole cular cloud. 

iMaddalena et al.l l| 19861 ) concluded that, in spite of the 
different velocities, Orion East is probably located at the 
same distance as Orion B, regarding its cometary shape and 
bright rim, pointing t oward Barnard's Lo op. This conclu- 
sion was questioned bv lKnude etaH (|2002h . who studied the 
distance to interstellar extinction features by using spectral 
types of the Michigan Catalog and Tycho-2 photometry, 
and found an absorbing cloud located at some 160 pc from 
the Sun toward the line of sight of L 1622, suggesting tha t 
this cloud may be a foreground object. IWilson et alJ l|2005h . 
based on the parallaxes of three Hipparcos stars, find that 
L 1622 is situated as close as 120 pc to the Sun. Nearby star 
forming clouds are of special interest because they allow us 
to study fine structural details and detect the lowest lumi- 
nosity objects. 

Motivated by these findings we started a spectroscopic 
and photometric study of the pre-main sequence stars and 
candidates associated with L 1622. Our aim is to assess the 
probable star forming history of the cloud, and see if the 
measured properties of these objects support the smaller or 
the larger distance. 

In order to have a homogeneous data set on the young 
low mass stars born in L 1622 we obtained moderate resolu- 
tion optical spectra and VRcIc photometry of all pre-main 
sequence stars and candidates projected on L 1622. Spitzer 
archive data were used for studying the spectral energy dis- 
tributions (SEDs) of the young stars. In order to find the 
physical properties of the cloud, compare the positions of 
the stars and the molecular g well as estimate the ef- 

ficiency of star formation we analysed 12 CO and 13 CO ob- 
servations of the cloud available in the data archive of the 
NANTEN radio telescope at Nagoya University. Our obser- 
vations and results are described in Sect. 2, and discussed 
in Sect. 3. Section 4 gives a short summary of our findings. 



2 OBSERVATIONS AND RESULTS 
2.1 Spectroscopy 

2.1.1 FAST spectra of PMS candidates 

Low-resolution optical spectra of the 16 Ha emission star s 
of the field of L 1622, listed by lOgura fc Hasegawal (|l983h . 
were acquired with FAST, a high throughput slit spectro- 
graph mounted at the Fred L. Whipple Observatory 1.5- 
m te lescope on Mount Hopkins, Arizona (|Fabricant et al.l 
1998). In addition to the Ha emission stars we observed 
HD 288313 too. We observed both component s of the visual 
binary L 1622-10 (|Reipurth fc Zinneckerll993h . separated by 
some 4.4 arcsec. 

We used a 300 g mm -1 grating blazed at 4750 A, a 
3-arcsec slit, and a thinned Loral 512 x 2688 CCD. These 
spectra cover 3800-7500 A at a resolution of ~ 6 A. We 
reduced and analysed the spectra in IRAfQ. 

After trimming the CCD frames at each end of the 
slit, we corrected for the bias level and flat-field for each 
frame, applied an illumination correction, and derived a 
full-wavelength solution from calibration lamps acquired 
immediately after each exposure. The wavelength solu- 
tion for each frame has a probable error of ±0.5-1.0 A. 
To construct final 1-D spectra, we extracted object and 
sky spectra using the optimal extraction algorithm within 
APEXTRACT. Most spectra have moderate signal-to-noise, 
S/N ~ 30 per pixel. The spectral classification was per- 
formed by visually comparing the spectra of MK standards 
with the s pectra of our target s follow ing the procedures de- 
scribed in ljaschek fc Jas chck (198?]). We estimate an accu- 
racy ±1 subclass for our spectral classification. For the LkHa 
stars our deriv e d spe ctral types agree with those found by 
ICohen fc Kuhil (|l979l ). The FAST spectra, normalized with 
the value of the continuum at 6560 A, are shown in Fig. [TJ 

2.1.2 CAFOS spectra of PMS candidates 

Part of the target stars have been observed with the CAFOS 
instrument on the 2.2-m telescope of Calar Alto Observa- 
tory in February 2007. Using the grism R-100, the observed 
part of the spectrum covered the wavelength interval 5800- 
9000 A. The spectral resolution of CAFOS observation, us- 
ing a 1.5-arcsec slit, was A/AA » 1000 at A = 8500 A. We 
reduced and analysed the spectra using standard IRAF rou- 
tines. The wavelength range of spectra was suitable for de- 
termining several flux ratios defined as too ls for spectral clas- 
sificat ion of late K and M type stars by Kirkpatrick et all 
i|l99ll ) ( A, B, C, B/A. B/C). TMartm fc .Kuril (Il996l) (h. 

and iPreibisch. Guenther. Zinneckerl (1200 ll ) (Tl, T2). 
We measured these spectral features on the spectra of our 
stars, and calibrated them against the spectral type and lu- 
minosity class by measuring them in a series of standard 
stars published bv lLe Borgne et al.1 (|2003l ). The accuracy of 
the spectral classification, estimated from the range of spec- 
tral types obtained from different flux ratios, is ±1 subclass 
(for further details of spectral classification see iKun et all 

1 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for the Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. |http:/ /iraf.noao.edu/| 



2004). The spectral types, obtained from the CAFOS obser- 
vations, are consistent with those obtained from the FAST 
spectra. 

We accepted an object as a classical T Tauri 
star if it fulfilled the cri terion established by 
iBarrado v Navascues fc Martini [|2003h . i.e. EW(Ra) 
exceeded the saturation limit of the chromospheric ac- 
tivity, depending on the spectral type. The Ha emission 
of such stars is thought to originate from the gas ac- 
creted by the star from the disc. The spectra of L 1622-1, 
L 1622-2, L 1622-7, L 1622-13, L 1622-14, and L 1622-16 
have shown no pre-main sequence signature. We no- 
ticed, however, that the spectrum of a faint M-type star, 
2MASS J 05540201+0140558, which appeared in the slit 
next to the CAFOS spectrum of L 1622-6, has shown 
W(Ha) « 15A. The critical EW(Ua) for its spec tral type 
M3, according to IBarrado v Navascues fc Martini , is about 

11 A, suggesting that 2MASS J 05540201+0140558 may be 
a CTTS. The equivalent width of the Ha emission line of 
M-type stars, however, may be slightly overestimated at 
the moderate resolution of our observations because, due 
to the TiO band and 65 69 A, the position of t he continuum 
is somewhat uncertain (|White fc B asri 2003). Taking into 
account this effect we regarded this star as a candidate 
pre-main sequence object and designated it as L 1622-6N. 
HD 288313 is a weak-line T Tauri star, whose pre-main 
seq uence nature is show n by the strong Li I absorption 
(cf. iMekkaden et~al]|2007h . Its spectrum shows a weak Ha 
emission, as the photospheric Ha absorption, expected for 
a K2 type star, is nearly filled in with emission (see the last 
panels of Figs. [T]and[2]). According to the £W(Ha) values, 

12 of 14 PMS stars covered by our spectroscopic survey are 
accreting objects. 

Table [T] shows the final list of the targets which have 
shown spectra characteristic of pre-main sequence stars. In 
addition to the derived spectral types we present the equiv- 
alent widths of the Ha and Li I lines in A, as well as indi- 
cate the additional emission lines observed in the spectra. 
The uncertainties given in parentheses have been derived 
from the repeatability of the measurements. We used the 
deblend command of the splot task of IRAF to separate 
the Li I line from the neighbouring absorption or emission 
features (CalA6718, [Sll] A 6717). The CAFOS spectra are 
shown in Fig. [2] The strong Ha emission lines show that all 
these objects are classical T Tauri stars. L 1622-3 exhibits an 
especially rich emission spectrum, characteristic of strongly 
accreting or eruptive (EXor type) stars. The spectral types 
of the st ars were converted int o effec tive temperatures T e a 
following iKenvon fc Hartmannl l|l995l ). 

2.1.3 Spectrum of VDB 63 

In addition to the candidate PMS stars we observed the 
spectrum of HD 288309, the illuminating star of the reflec- 
tion nebula VDB 63, located at the eastern edge of L 1622 
(see Fig. [6)1, on the 4900-7800 A region, using CAFOS with 
the grism G-100. The spectral type and photometric data 
of this star may constrain the distance of L 1622. We de- 
rived a spectral type K0 HI. The abso l ute m agnitudes of K0 
giants, tabulated bv lWainscoat et~aL 

I (|1992| ). together with 



the photometric data of HD 288309 found in the NOMAE0 
catalogue, result in a most probable distance D ~ 320 pc for 
VDB 63. The large di spersion of the absolu te magnitudes 
of K giants, tjyi = 0.7 (|Wainscoat et al.lll992T ). lends a high 
uncertainty to this value. Taking into account the 1-cr range 
of absolute magnitudes, the corresponding distance interval 
is 230 pc < D < 440 pc. 

2.2 VRcIc photometry 

Photometric observations of the young stars in the V, Rc 
and Ic bands were undertaken on 15 January and 14 Febru- 
ary 2007 using the 1-m RCC-telescope of Konkoly Obser- 
vatory. We used a Princeton Instruments VersArray:1300B 
camera, that utilizes a back-illuminated, 1300x1340 pixel 
Roper Scientific CCD. The pixel size is 20 fim, correspond- 
ing to 0.31 arcsec on the sky. Integration times were between 
5 s and 300 s. In order to calibrate the observed fluxes, we 
bracketed the o bservation of ea ch target star by observations 
of four nearby lLandoltl (1 19831 ) standard stars. We reduced 
the images in IRAF. After bias subtraction and Hatfield cor- 
rection psf photometry was performed using the phot and 
allstar tasks of the DAOPHOT package. 

The results of the photometry are presented in Table [2] 
The photometric errors, given in parentheses, are quadratic 
sums of the formal errors of the instrumental magnitudes 
and those of the coefficients of the transformation equations. 

2.2.1 Interstellar extinction 

The unredde ned colour indices of the sp ectral types were 
adopted from lKenvon fc Hartmannl (1 19951 ). The interstellar 
extinction suffered by the program stars can be derived from 
their colour excesses. The colour excesses of pre-main se- 
quence stars, however, reflect not only the interstellar red- 
dening, but are also influenced by several other effects, such 
as the light, scattered from the disc and envelope surround- 
ing the star, as well as thermal emission from the disc and 
accretion shocks on the stellar surface. It is usually assumed 
that most of the E n n -i r , colour excess results from interste l- 
lar reddening (e.g lHartigan et al. I ll994l ; iMever et al.lll997l ), 
and thus it is a reliable estimate of the extinction usin g the 
relationship Ay 1 = 4.76 x E Rc -i c l|Cohen et al.lll98ll ). 

During our data analysis, however, we realized that 
the Er c -i c colour excesses of LkHa336c, L1622-10A, 
and L1622-10B are probably influenced by scattered light: 
Er c -i c apparently underestimates the extinction of these 
stars: the luminosities, derived from the Ic magnitudes, cor- 
rected for the extinction based on Er c -i c , are too low for 
classical T Tauri stars, located in the distance interval 120- 
400 pc from the Sun. The presence of scattered light is sup- 
ported by our photometric images, which show the profiles 
of these stars slightly broader than the average stellar pro- 
file of the image. We therefore determined the extinction 
using the Ic — J colour excess as well, using the 2MASS 
|Cutri et al.ll2003T ) data, and compared the extinctions ob- 
tained from both methods. We found that in most cases 
both colour excesses result in virtually the same Av, sug- 
gesting that both the influence of scattered light in the Rc 
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Table 1. Results of spectroscopy 



Star 


Sp.T. 


-EW (Ha 

(A) 


) EW{Ul) 

(A) 


Other emission lines 




K2 


0.80 (0.1) 






J_j _L<i — _ O 


MO 


zyj . u . \j } 




Rnlmpr Gfvi^G f^n TT TT \C- C^a TT TR tviT^lnt 
l_>dlllllJI V_jd 11 1± Oc. I\ 5 vd 11 11\. L11JJ1UL, 

HeIAA5873, 6678, Nal, numerous Fell lines 


LkHa 334 


K4 


-8.5(0.5) 


0.65 (0.03) 


H/3, H7; He 1 A 5873, Can IR triplet 


LkHa 335 


K4 


-52.5(2.0) 




Balmer series, He 1 AA 5873, 6678; Ca II IR triplet 


L 1622-6 


Ml 


-20.5(1.0) 


0.51 (0.05) 


H/3 


L 1622-6N 


M3 


-14: 






L 1622-8 


M2.5 


-103.0(1.0) 




Balmer series; Call H & K; He 1 A 5873 


LkHa 336 


K7 


-26.0(1.0) 


0.29 (0.05) 


Balmer series, Can H & K 


LkHa 336c 


Ml 


-54.8 (2.0) 


0.53 (0.05) 


H0, H5 


L 1622-10 A 


K7 


-41.0(3.0) 




[Oi] AA 6300, 6363; [Nil] A 6583, [S 11] AA 6716, 6731 


L 1622-10 B 


K5 


-79.9(3.0) 




H/3 


L 1622-11 


M2 


-65.1 (1.0) 


0.81 (0.02) 


Balmer series, Can H & K 


LkHa 337 


K7 


-50.5(1.0) 


0.45 (0.05) 


Balmer series, Call H & K; He 1 AA 5873, 6678; [Oi] A 6300 


L 1622-15 


MO 


-72.2(2.0) 




Balmer series, Call H & K; MglA5172; [0 1] AA 6300, 6363 



He 1 AA 5873, 6678; [S 11] AA 6716, 6731; Fen 



band and the J-band emission from the disc inner rim are 
negligible. Exceptions are the above stars. Probably scat- 
tered light makes these stars bluer in the Rc — Ic, and, 
at the same time, their Ei c -j colour excess is slightly in- 
creased by the po ssible contribution of the disc emission to 
the J magnitude l|Cieza et al.ll2005T ). Therefore we adopt as 
a best estimate of Ay the average of A v l = 4.76 x Er c -i c 
and Ay = 3.23 x Ei c -j for each star. The Ay values, as 
well as the 2MASS designations of the target stars, formed 
from their equatorial coordinates, are also given in Ta- 
ble [5] For determining the extinction in the other photo- 
metric bands, we used the relati ons An r , = 0.78 A y and 
A Ic = 0.59 Ay for the optical (| Cohen et all Il98lf ), and 
A j = 0.26 Ay, A H = 0-15 Ay, A Kb = 0.10 Ay, presented 
bv lNielbock fc Chinil (|2005h for the 2MASS-bands. 

Bolometric luminosities were derived b y applying the 
bolom etric corrections BCi c , tabulated by IHartigan et alj 
(|1994[ ) to the dereddened 7c magnitudes, based on the as- 
sumption that the total emission of the s tars in the Ic band 
originates from the photosphere (e.g. iMever et al. I 1 19971 ; 
ICieza et al.ll2005h . We determined Z/ bo i/LQ for various dis- 
tances found in the literature. 



2.3 Infrared archive data 

2.3.1 IRAC observations 

We measured the IRAC fluxes of the program stars using 
images from Spitzer Space Observatory archive (Prog ID 43: 
An IRAC Survey of the L 1630 and L 1641 (Orion) Molec- 
ular Clouds - Fazio et al.). We created the IRAC mosaics 
from Basic Calibrated Data (BCD) frames in each chan- 
nel using a custom IDL script which also converted the 
MJysr -1 units to DNs -1 using conversion factors 0.1088, 
0.1388, 0.5952, and 0.2021 for channels 1, 2, 3, and 4, re- 
spectively, and performed aperture photometry on all im- 
ages with Phot Vis version 1.10 , which is an IDL-GUI based 
photometry visualization tool (|Gutermuth et al] [2004l . The 
radii of the source aperture and the inner and outer sky 
annuli were at 2.4 arcsec, 2.4 arcsec, and 7.2 arcsec, respec- 



tively. We calculated the standard IRAC magnitudes us- 
ing —2.5 log(flux£(jv/s)+ z P, where the zero-point term, zp, is 
21.9929, 21.2583, 19.0747, and 19.4372 for channels 1-4, re- 
spectively. The zero point includes the zero point of the flux 
to magnitude transformation based on the large aperture 
measurements of a standard star (19.66, 18.94, 16.88, and 
17.39 mag for channels 1-4); and the aperture corrections to 
account for the difference between the aperture sizes used for 
the standard star and for L 1622 photometry ( —0.21, —0.23, 
—0.35 , and —0.50 mag for channels 1-4). See iReach et al.1 
(2005) for the detailed description of the calibration of the 
IRAC data. We converted the magnitude into flux units us- 
ing the equation / = c x io~°' 4 * m where / is the flux m is 
the corresponding magnitude and c is the zero point in Jy 
(280.9, 179.7, 115.0, 64.13 for channels 1, 2, 3, and 4 respec- 
tively). The resulting IRAC magnitudes are listed in Table[3] 
L 1622-3, L 1622-15, and LkHa 334 are located outside the 
field of view of the IRAC observations. 



2.3.2 MIPS data 

We also present observations of L 1622 at 24 and 70 /im, ob- 
tained with the Multiband Imaging Photometer for Spitzer 
(MIPS) as part of program PID 43 (AOR 16214784, Nov 4 
2005). The data were downloaded from the Spitzer Science 
Center Archive with Leopard. The cloud was imaged using 
a single scan map with 4 legs at full array cross-scan offset, 
resulting in half array gaps at 70 ^m. The scanlegs were 
half degree in length. Observations were carried out using 
medium scan mode, resulting in a total effective exposure 
time of 40 s/pixel at 24, and 20 s/pixel at 70 /im. 

All data were processed usin g the MIPS instrum ent 
team Data Analysis Tool (DA T IGordon etafl 120051) as 
descri bed by lEngelbracht et al"] (|2007 ) and IGordon et al.l 
|2005l ). Care was taken to minimize instrumental artifacts. 
We used the MIPS calibration star HD 173398 as a PSF 
standard for the 24 (im photometry, with the final PSF con- 
structed from 72 individual AORs, kindly provided to us by 
C. Engelbracht. 

Standard IRAF tasks phot and allstar of the DAOPHOT 
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Figure 1. FAST spectra of the pre-main sequence stars associated with L 1622. The separate panel at the bottom shows the part of the 
spectrum of HD 288313 around the Ho and Li I. 



Table 2. VRqIc magnitudes and derived visual extinctions of the pre-main sequence stars of 
L1622 



Star 


2MASSJ 




V 






Rc 




Ic 


A v 


HD 288313 


05540300+0140218 


10 


014(0 


02) 


9.263 (0.02) 


8.828(0.02) 





93 


L 1622-3 


05533298+0217200 


16 


365 (0 


02) 


15 


282 (0.02) 


14 


171 (0.02) 


1 


33 


LkHo 334 


05534090+0138140 


13 


122 (0 


02) 


12 


365 (0.02) 


11 


617(0.02) 


1 


48 


LkHo 335 


05535869+0144094 


14 


201 (0 


02) 


13 


164 (0.02) 


12 


276 (0.02) 


2 


57 


L 1622-6 


05540189+0140272 


15 


862 (0 


02) 


14 


464 (0.02) 


12 


635 (0.02) 


2 


45 


L 1622-6 N 


05540201+0140558 


18 


421 (0 


06) 


16 


866 (0.04) 


14 


882 (0.03) 


3 


24 


L 1622-8 


05540718+0139007 


17 


424 (0 


05) 


15 


190 (0.04) 


13 


584 (0.03) 


1 


88 


LkHo 336 


05542008+0142565 


14 


229 (0 


04) 


13 


188 (0.04) 


12 


278 (0.03) 


1 


70 


LkHo 336c 


05542031+0142518 


18 


160 (0 


08) 


16 


707 (0.08) 


15 


416(0.06) 


3 


67 


L 1622-10 A 


05542669+0152154 


17 


899 (0 


08) 


16 


467 (0.08) 


14 


857(0.06) 


4 


83 


L 1622-10 B 


05542696+0152172 


18 


162 (0 


08) 


16 


768 (0.06) 


15 


457(0.05) 


6 


05 


L 1622-11 


05543040+0150593 


17 


596 (0 


05) 


16 


174 (0.03) 


14 


337 (0.03) 


1 


80 


LkHo 337 


05543793+0129513 


14 


406 (0 


02) 


13 


589 (0.02) 


12 


507 (0.02) 





72 


L 1622-15 


05571611+0229074 


15 


650 (0 


02) 


14 


974 (0.02) 


13 


950 (0.02) 


1 


66 
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Figure 2. CAFOS spectra of pre-main sequence stars associated with L 1622. The separate panel at the bottom shows the part of the 
spectrum of HD 288313 around the Ha and Li I. 



Table 3. IRAC and MIPS magnitudes of the pre-main sequence stars of L 1622 



Star 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 


LkHa 334 










3.95 (0.01) 


LkHa 335 


7.290 (0.001) 


6.783 (0.002) 


6.287(0.001) 


5.289 (0.001) 


2.52(0.01) 


L 1622-6 


9.316(0.004) 


9.147 (0.003) 


8.967 (0.005) 


8.492 (0.006) 


3.88 (0.01) 


L 1622-6 N 


10.753 (0.003) 


10.335 (0.003) 


9.971 (0.005) 


9.226 (0.005) 


6.78(0.01) 


L 1622-8 


9.827(0.002) 


9.470 (0.002) 


8.934 (0.003) 


8.069 (0.002) 


4.96(0.01) 


LkHa 336 


8.374 (0.003) 


7.969 (0.002) 


7.542 (0.002) 


6.541 (0.004) 


3.52 (0.02) 


LkHa 336c 


9.612(0.005) 


9.196 (0.005) 


8.632(0.003) 


7.644 (0.003) 


2.95 (0.02) 


L 1622-10 A 


8.838 (0.002) 


8.358 (0.002) 


7.817(0.004) 


7.068 (0.002) 


3.11 (0.02) 


L 1622-10 B 


9.620 (0.010) 


9.160,(0.010) 


8.570 (0.020) 


7.700 (0.010) 


3.00 (0.02) 


L 1622-11 


10.618(0.002) 


10.289 (0.003) 


9.926 (0.005) 


9.201 (0.004) 


6.36(0.01) 


LkHa 337 


9.679 (0.002) 


8.771 (0.003) 


8.577 (0.003) 


7.548 (0.002) 


4.33 (0.01) 


IRAS 05522+0146 


12.491 (0.005) 


11.190 (0.004) 


10.056 (0.005) 


9.095 (0.003) 


3.69 (0.01) 



package were used to carry out the 24 fim PSF and 
70 /im aperture photometry. Since the MIPS psf is 
very stable and the nebulosity effects were negligible at 
24 /xm, we used a large PSF radius of 112 arcsec, with 
fitting radius of 5.7 arcsec. The instrumental number 
flux counts were converted to fluxes with the conversion 
1.068x10" 3 mjy arcsec" 2 MIPS.UNIT" 1 as determined by 
the MIPS team (|Engelbracht et al.l [2007h . We then trans- 
lated these to [24] magnitudes taking 7.17 Jy for the [24] 
magnitude zero point. 

The large PSF radius used ensured us that the aperture 
correction was smaller than 0.06 mag at [24] (the aperture 
corre ction for 35-arcsec aperture radius; [Engclbrach t et ail 
|2007| ). The correction factor for the PSF used has been de- 
termined by the MIPS team during previous work (Caspar 
et al. 2008, in prep.) to be only 0.0368 magnitudes (3.4%). 



We summarize our 24 photometry results in Tabled L 1622- 
3 and L 1622-15 are located outside the field of view of the 
MIPS observations. Three stars could be identified in the 
70 /im images, namely LkHa 335, LkHa 336c, and L 1622- 
10A. We performed aperture photometry for these stars us- 
ing a 18 arcsec aperture radius. The sky contamination was 
estimated in an annulus outside the aperture with inner and 
outer radius of 18 arcsec and 36 arcsec respectively. We con- 
verted the MIPS unit to Jy using a calibration factor of 
16.52 mjy/square arcsec resulting in fluxes of 529±341 mjy, 
927±404 mjy, and 504±302 mjy for LkHa 335, LkHal- 
pha336c, and [UTF-8?]L 1622-10A.respectively. No useful 
stellar data could be extracted from the 160 /im image. 




Figure 3. Left: The 12 CO spectrum observed at 
(1, 6)=(204?8, -11?8); Right: The average 12 CO spectrum 
of L 1622. 



2.4 12 CO and 13 CO maps of L 1622 

Observations were carried out in the 12 CO(J = 1-0) emis- 
sion toward the Orion B region with the 4-m millimeter- 
wave telescope, NANTEN, at the Las Campanas Observa- 
tory, from November 2001 to March 2002. The half-power 
beam-width (HPBW) was 2.6 arcmin at 115 GHz. The front- 
end was a superconducting mixer followed by HEMT ampli- 
fiers. D etails of this mixer receiver are given bv lOgawa et al.l 
(1990). The backend was an acousto-optical spectrometer 
with a velocity resolution of 0.1 kms -1 . The telescope point- 
ing was measured to be accurate within 20 arcsec by observ- 
ing planets. The spectral intensities were calibrated by em- 
ploying th£^tandardroorn-^temperature chopper-wheel tech- 
nique (|Kutner fc Ulicbl [l981 V The absolute intensity cali- 
bration was made by observing Orion-KL [R.A.(B1950) = 
5 h 32 m 47.5 s , Decl. (B1950) = -5°24'21"], by assuming its 
absolute temperature, ~ 65 K. We also observed it ev- 
ery 2 hr to confirm the stability of the whole system. Most 
of the observations were made by the frequency switching 
technique with a frequency interval of 15 MHz, while the 
position switching technique was also used when the atmo- 
spheric frequency is close to the celestial 12 CO frequency. 

These observations were made as a part of the large 
scale 12 CO survey of the Orion-Eridanus region which con- 
sists of ~ 100,000 observing points. In the Orion region, 
b > —30°, data were obtained at a grid spacing 4 arcmin. 
The r.m.s. noise fluctuations were ~ 0.53 K per channel for 
a typical integration time of ~ 8 s. 

The 13 CO (J = 1-0) observations were made at 4- 
arcmin grid spacing toward the same region as a part of 
large scale 13 CO survey of the Orion region in January and 
February 2002. The half-power beam-width (HPBW) was 
2.7 arcmin at 110 GHz, and the r.m.s. noise fluctuations 
were ~ 0.16K per channel. The frequency-switching tech- 
nique with a frequency interval of 15 MHz was also used. 
We observed Orion-KL every 2 hr to check the stability and 
establish the absolute intensity scale. We assumed its abso- 
lute temperature Tr ~ 10 K. 

The left panel of Fig. [3] shows the 12 CO spectrum ob- 
served at (J,&)=(204°8, -11?8). The self-absorbed feature at 
8-13 kms -1 corresponds to Orion B, whereas the peak at 
lkms -1 represents L 1622. In the right panel the average 
of all 12 CO spectra within the 3-ct boundary of the 12 CO 
cloud, corresponding to L 1622, is plotted. The peak veloc- 
ity of this spectrum is +1.17 kms -1 , and the line width at 
half maximum is Av — 1.66 kms - . 
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Figure 4. Large-scale distribution of the molecular gas in the 
region of Orion B. Left: 12 CO emission integrated over the veloc- 
ity interval 5 kms" 1 < ijlsr < 15 kms" 1 for the Orion B main 
cloud. Right: — 5kms -1 < i>lsr < +5kms -1 map, containing 
the L 1622 molecular cloud. St ar symbols ind i cate t he 0-B3 type 
members of Orion OB lb from iBrown et ah 

I dl994l') , to show the 
relation between molecular clouds and Ori OBlb. 



The large-scale distribution of the molecular gas is dis- 
played in Fig. 2] The left panel shows the distribution of 
the 12 CO emission integrated over the velocity interval of 
(5,15) kms -1 , and the right panel shows the same distri- 
bution over the velocity interval ( — 5, +5.0) kms - . Using 
the conversi on formula JV(H 2 ) / W{CO) = 1.8xl0 20 cm -2 

KB), 



/ Kkms -1 (|Wilson et al.ll2005l ), our CO data allow us to 
estimate the mass of the cloud. Summing up the obtained 
column densities within the 3-a contour (1.5 Kkms -1 ) of 
the cloud, and adopting a mean molecular weight of 2.333, 
the result is M c ioud ~ 1100 Mq, if we assume a distance 
of 400 pc. At a distance of 120 pc the measurements re- 
sult in a mass of 100 Mq . The virial mass of the cloud, 
Mvm = 190 x R x Av 2 « 1500 Mq at 400 pc, where R is 
the average radius derived from the 3-cr area of the cloud. 
The mass of the 13 CO emitting gas is M L te( 13 CO) = 



750 Mq, and 



M V ir( 13 CO) = 



770 Mq. Thus the ratio 



M( 13 C O) / M( 12 CO) « 70%, a typical value for nearby dark 



cloud s ( de Geus. Bronfman. Thaddeusl Il990l ; iMizuno et al.l 
1 19951 . Il998h . Mvtr /Mt.tr ~ 1 for both the 12 CO and the 
13 CO clouds. 

The column density of the molecular hydrogen to- 
ward the 13 CO peak at (l,b) = (204.8,-11.8) is N(H 2 ) = 
9.3 x 10 21 cm -2 . The typical column densities of the isolated 
small 13 CO clouds around the main cloud are about N(H 2 ) = 
1.9 x 10 21 cm -2 . Figure [5] shows the 12 CO and 13 CO maps 
of L 1622, with the positions of the associated PMS stars 
overplotted. 



3 DISCUSSION 

3.1 Description of the region 

L 1622 harbours two IRAS point sources whose flux ratios 
are different from those of ordinary stars. IRAS 05517+0151 
coincides with a visible star. Both the flux ratio F12/-F25 
and the near-infrared properties of its associated 2MASS 
source, 2MASS J 05542277+0152039, are indicative of an 
evolved star. The optically invisible IRAS 05522+0146 
is probably an embedded young stellar object. Within 
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Figure 5. Left: The new 12 CO map of L 1622, with the posi- 
tions of PMS stars studied in this paper, and a grid of equa- 
torial (J2000) coordinates overplotted. The lowest contour is at 
1.5 K kms -1 , and the increment is 5Kkms -1 . Left: 13 CO map 
of L 1622. The lowest contour and the steps are IK kms -1 and 
2Kkms -1 , respectively. The emission is integrated over the ve- 
locity interval — 5kms _1 < V+SR < +5kms _1 . in both maps. 

the uncertainty ellipse it coincides in position with 
2MASS J 05545377+0147098. This source was not detected 
at 12 fim by IRAS, but the high upper limit of the 12 )im 
flux, 0.34 Jy, is compatible with the assumption that the 
IRAS and 2MASS fluxes originate from the same protostel- 
lar source. No infr ared source can be found at the position 
of HH 122 VLA1 (|Rodriguez fc Reipurth|[l994T ). 

Figure [6] shows the surface distribution of the 12 CO 
emission and star formation signposts, overplotted on the 
DSS 2 red image of L 1622. The value of the lowest 12 CO con- 
tour plotted is 1.5 K kms" 1 , and the increment is 5 K kms" 1 . 
Crosses indicate the pre-main sequence stars, and dotted e l- 
lipses show the optically selected cores l|Lee fc Mv crs 1999). 
Positions of HH 122, HH 122 VLA 1, and IRAS 05522+0146, 
as well as the reflection nebulae VDB 62 and VDB 63 are 
also indicated. 

The 12 CO distribution shows a cometary shaped cloud 
with a bright rimmed head region on the southwestern side. 
This shape suggests that star formation might have been 
triggered by a shock arriving from the southwest, where 
the luminous stars of Ori OBI can be found. Four pre- 
main sequence stars, namely L 1622-3, LkHa 334, LkHa 337, 
and L 1622-15 are projected outside the observed molecular 
cloud, thus their relation to L 1622 is uncertain. In particu- 
lar, L 1622-3 and L 1622-15, located well outside the lowest 
contour of the tail of the cloud, may be unrelated to L 1622. 
LkHa 334 and LkHa 337 are projected near the head of the 
cloud, on the southwestern side, and thus might have been 
born in L 1622 at an earlier epoch when the cloud was larger. 
On the other hand, these isolated stars might have been 
formed in small clouds around the main large cloud, as dis- 
cussed in iMizuno et alj (Il998l ) for the Chamaeleon region. 
Ten PMS stars are projected within the head of the cloud. 
Four of them, HD 288313, L 1622-6, L 1622-8 and L 1622- 
6N form a compact group near the bright rim. Three stars, 
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Figure 6. Distribution of the 12 CO emission and the young stel- 
lar objects, overplotted on the DSS2 red image of the L 1622 re- 
gion. The value of the lowest contour plotted is 1.5 K kms -1 , 
and the increment is 5Kkms -1 . Crosses indicate the pre-main 
seque nce stars, and dotte d ellipses show the optically selected 
cores l lLee fc Mversl Il999h . Position of HH 122, HH 122 VLA1, 
and IRAS 05522+0146, as well as the reflection nebulae VDB 62 
and VDB 63 are also indicated. The beam size of the 12 CO ob- 
servations is shown by the circle at the lower left corner. 

L 1622-10 A, L 1622-10 B, and L 1622-11, are projected on a 
dark spot within the cloud, associated with the 13 CO peak. 
Both candidate embedded protostars, HH 122 VLA 1 and 
IRAS 05522+0146, are associated with the core L 1622 B. 

3.2 Positions of the pre-main sequence stars in 
the HRD 

The Hertzsprung-Russell diagrams (HRD) of the pre-main 
sequence stars of L 1622, for four different distances, are 
displayed in Fig. [7] Error bars of log T c h are derived from 
the accuracy of the spectral classification. The error of the 
luminosity comes from the quadratic sum of 81c, SAi c , 
and 5BCi c . Evolutionary tracks and isochrones, as well 
as the position of the birthline and zero-age main-sequence 
(|Palla fc Stahlerlll999h are also shown. 

Most of the programme stars are distributed over a 
fairly narrow age range. Exceptions are HD 288313 and 
LI 622-6, which are apparently far more luminous than 
the other stars. Sinc e HD 288313 is a multiple system 
ijReipurth et al-lbOOSal ). the observed luminosity has to be 
distributed on its components. L1622-6 is either a fore- 
ground object, or, more probably, a multiple system as well. 
The distribution of the stars with respect the isochones 
shows that the extreme emission line star L 1622-3 is appar- 
ently the oldest member of the sample, in contrast with the 
observations that PMS activities decrease with stellar age. 
It is unlikely that the relatively low luminosity of this star is 
due to its larger age. L 1622-3 may either be a higher-mass 
counterpart of the underluminous, highly accreting very low 
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Figure 7. Positions of the observed stars in the HRD assuming that the distance of the cloud is (a) 120 pc, (b) 160 pc, (c) 320 pc, (d) 
400 pc. Thin solid lines indicate the isochro nes as labelled, and dot ted lines show the evolutionary tracks for the masses (Mq) indicated 
at the lower end of the tracks according to iPalla fc Stahlerl |l99§) model. The dash-dotted line corresponds to the birthline and thick 
solid line indicates the zero age main sequence. 



mass stars LS-RCr A 1. Par-Lup3-4 dComeron et alj 120031 ) 
and KPNO Tau 12 jLuhman et all 120031 ') . or may be some- 
what more distant than L 1622. Its position on the sky with 
respect to the cloud suggests that L 1622-3 might have been 
born in another small cloud, projected near L 1622. Ra- 
dial velocity measurements may clarify its connection with 

L1622. 

At a distance of 120 pc (panel fal lWilson et alj|200d ). 
most stars are close to the ZAMS, while L1622-3 and L1622- 
15 are on the ZAMS. The average age of the stars, suggested 
by this distance, is not compatible with their CTTS nature 
and triggering of star formation in L 1622 from the older 
subgroups of Ori OBI. The same is true for 160 pc, shown 
in panel (b): the mean age of the stars would be about 20 
million years at this distance. At 320 pc (panel (c)), sug- 
gested by the photometric properties of VDB 63, the mean 
age of the pre-main sequence stars is about 3 million years, 
whereas, if L 1622 is situated at 400 pc from us, the age 
of its YSO population is around 1 million years (panel (d) 
of Fig. |7J). We propose to accept 400 pc for the distance of 
L 1622, and thus the age of < 1 million year for the aver- 
age age of its PMS star population. This age is compatible 
with that suggested by the spectroscopic properties of the 
stars, i.e. the high pr oportion of accretors (85% ) among the 
Ha emission stars (cf. lJavawardhana et alll2006l ). as well the 
high propor tion of the protostars among the YSO popula- 
tion (10/28, iReipurth et al.ir2008bl ). 



3.3 Spectral energy distributions of the target 
stars 

The shape of the SED, reflecting the structure of the circum- 
stellar matter of the YSOs, is a useful indicator of evolution- 
ary s tate (e.g lLada fc Wilkind Il984l ; lAdams. Lada, fc Shul 
1 19871 ). The optical, near- and mid-infrared parts of the 
spectral energy distributions of the programme stars were 
constructed using the V RcIcJ H K s magnitudes and IRAC 
and MIPS fluxes, each corrected for the interstellar extinc- 
tion. The fl uxes correspo nding to zero magnitude were ob- 
tained from iGlassI (|l999h for the VRcIc bands, and from 
the 2MASS All Sky Data release web document for the 
JHK a bands. IRAS fl uxes are also plotted when available 
l|Weaver fc Jones! Il992l ). The resulting SEDs are shown in 
Fig-El Photospheric SEDs have been drawn by dashed lines. 
Their values were determined from the dereddened Ic mag- 
nitudes, and from the colour indices corresponding to the 
spectral types. The right bottom panel of Fig. [8] shows the 
fluxes of IRAS 05522+0146, with the assumption that it is 
identical with 2MASS J 05545377+0147098. In this case the 
fluxes are not corrected for the interstellar extinction. The 
tip of the downward arrow indicates the upper limit of its 
12- micron flux. 

The SEDs exhibit near-infrared (H , K s ) colour excess 
of LkHa335, LkHa336, and L 1622-10 A, indicative of disc 
inner edges at the dust sublimation radius, characteristic of 
actively accreting stars, probably not older than 1 million 



3 http : //www. ipac.caltech.edu/2mass/releases/allsky/doc/ 
sec6_4a.html 
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Figure 8. Spectral energy distributions of the PMS stars and IRAS 05522+0146, based on our optical photometry (plusses), 2MASS 
(dots), Spitzer IRAC and MIPS (triangles), as well as IRAS (thin diamonds) data. The arrow in the last panel indicates the upper limit 
of the 12-micron IRAS flux of the source. 



years. The signature of the disc appears at longer wave- 
lengths in the SEDs of L 1622-6, L 1622-6N, L 1622-8, and 
L 16 22-11, suggesting m ore evolved discs around these stars 
fe.g lFurlan et al]|2006t ). We note that these latter stars, ex- 
cept L 1622-11 are members of the small compact group near 
the bright rim, whereas most of the stars with NIR excess 
are projected deeper inside the molecular cloud, suggesting 
star formation propagating into the cloud from the direction 
of the bright rim. We also note that the signature of the disc 
in the SED of L 1622-6N confirms the PMS nature of this 
star. 

In order to get a deeper insight into the circumstel- 
lar structures of the program stars, we fitted the observed 
SEDs with YSO model SEDs using the large grid of r adiation 
transfer models described bv lRobitaille et all (|2Q07tFl Since 
several combinations of central star-disc-infalling envelope 
can result in good fit with the observed SED, we selected 
those model SEDs from the best fits that reproduced the 
stellar parameter values (Av, T e g, M*) determined from our 
optical observations. Table [3] shows the main stellar and cir- 
cumstellar parameters resulted from the model fitting. The 
stellar parameters obtained from the optical observations 
are also shown for comparison. We have not fitted the SEDs 
of stars without Spitzer data. 

It can be seen that the best models for L 1622-6N, 
L 1622-8, L 1622-11, and LkHa 337 contain only a star and 
an accretion disc, whereas the SEDs of LkHa 335, L 1622- 
6, L1622-10A, LkHa 336, and LkHa 336c could not be fit- 

4 Available online at http://www.astro.wisc.edu/protostars 



ted without including an envelope. The shape of the SED 
of L 1622-6 ranks this object into the class of transitional 
discs, having a wide inner hole. Envelopes usually disappear 
before development of the inner disc hole, therefore we do 
not expect an envelope around objects whose SED shapes 
are similar to that of L 1622-6. The large inner disc hole may 
be a consequence of the binarity of this star, also suggested 
by its high luminosity in the HRD. The inner regions of the 
accretion disc of this object might have been evacuated by 
a close companion. 

In the case of LkHa 336c and L 1622-10A the fits sug- 
gest discs seen nearly edge on. Contribution of scattered 
light is significant in the short-wavelength part of the SEDs 
of these stars, as we have seen in Sect. 12.21 The SED of 
IRAS 05522+0146 shows this object to be a Class I source, 
with signs of strong water ice and silicate absorptions around 
3 and 10 /im, respectively. 

3.4 A possible scenario of star formation in L 1622 

The Hertzsprung-Russell diagrams presented in Fig. [7] sug- 
gest that L 1622 is not associated with the near side of the 
Orion-Eridanus Bubble. The extinction layers identified at 
120 and 160 pc, that represent the interaction region of the 
Local Bubble and Ori-Eri Bubble, are foreground to this 
cloud. The derived parameters of young stars are compatible 
with both the distance of Ori OBla (330 pc) and Ori OBlb 
(440 pc). The cloud is projected on Orion B, but is diverg- 
ing from this giant cloud toward us with a radial velocity of 
9 kms -1 . Thus the true distance of the cloud from the Sun 
is probably slightly smaller than 400 pc. This geometry sug- 



Table 4. Parameters of the programme stars, derived from SED model fitting and from our optic al data, assuming a dis tance of 400 pc, 
as well as the main properties of the circumstellar matter, according to the best fit models of iRobitaille et al.l l|2007f ). Stars without 
Spitzer data are not listed. 
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3405 


3525 


0.31 


0.35 


1.5 


;$ 1-° 






8.8E-3 


8.6E-8 


1.0 


0.04 


103 


32 


LkHa 336 


4120 


4060 


0.80 


0.80 


1.0 


< 1.0 


7.2E-9 


1.9E+3 


3.2E-2 


6.8E-7 


1.4 


0.80 


100 


76 


LkHa 336c 


3980 


3720 


0.70 


0.50 


0.2 


2.5 


9.4E-7 


4.0E+3 


4.1E-3 


2.4E-7 


1.0 


0.16 


56 


81 


L 1622-10 A 


4010 


4060 


0.70 


0.80 


0.4 


1.5 


4.0E-7 


1.2E+3 


2.8E-4 


5.7E-10 


1.0 


0.14 


115 


81 


L 1622-10 B 


4200 


4350 


1.01 


0.90 


0.3 


3.0 


9.3E-7 


1.6E+3 


6.0E-5 


2.0E-11 


1.0 


0.17 


984 


87 


L 1622-11 


3540 


3580 


0.37 


0.35 


1.9 


Si-o 






8.2E-5 


1.6E-11 


2.9 


0.11 


100 


32 


LkHa 337 


4030 


4060 


0.75 


0.80 


4.0 


1.5 






9.1E-3 


3.4E-8 


4.9 


0.30 


198 


69 



Column 1: Name of the star; Cols. 2—3: T e g derived from the SED model fit and spectral classification, respectively; Cols. 4-5: Mass of 
the star in M , derived from the SED model fit and HRD, respectively; Cols. 6-7: age of the star, derived from the SED model and 
HRD, respectively; Col. 8: accretion rate from the envelope in M yr - 1 ; Col. 9: outer radius of the envelope in AU; Col. 10: disc mass 
in M ; Col. 11: disc accretion rate in M yr - 1 ; Col. 12: inner radius of the disc as a multiple of the dust sublimation radius i? su b; 
Col. 13: same in AU; Col. 14: outer radius of the disc in AU; Col. 15: inclination (deg). 

Assuming an age difference of ~ 6 million years the average 
expansion velocity of the Hn zone or stellar wind bubble 
should have bee n about 6.8 kms -1 , a reasonable value fo r 
such effects (e.g lElmegreen fc Ladalll977l ; lElmegreenlll998u . 
Similar interaction between L 1622 and Ori OBla, assuming 
a distance of ~ 320 pc, is less likely, considering the orienta- 
tion of the bright rim of L 1622 and the projected position 
of Ori OBla. 



4 CONCLUSIONS 

We determined effective temperatures and luminosities of 
14 pre-main sequence stars and candidates in the region of 
the dark cloud L 1622 and found that these quantities are 
incompatible with the assumption that L 1622 lies at the 
near side of the Orion-Eridanus Bubble, but suggest that 
this cloud is located at the distance of Orion B. The derived 
ages of the stars at this distance are around 1 million years. 
The shapes of the SEDs of the target stars support this age 
estimate. Star formation in L 1622 was probably triggered 
by the luminous stars of the association subgroup Ori OBlb 
(Orion's Belt). 

We present new 12 CO and 13 CO maps of L 1622. Ac- 
cepting a distance of 400 pc we derived a mass of 1100 Mq 
for L 1622 from the 12 CO observations. The total mass of the 
YSOs identified so far, assuming a mean mass of 0.5 Mq 
for the Spitzer sources not included in our spectroscopic 
survey, is «20Mq, suggesting a star formation efficiency 
SFE w 1.8%. 

Several of our programme stars may be interesting tar- 
gets for more detailed studies. The derived luminosity of 
the CTTS L 1622-6 suggests that this object is a binary or 
multiple system. The large inner hole in the accretion disc 
of this star, suggested by the models fitted to its SED, sup- 
ports this hypothesis. The SEDs of L 1622-10A, L 1622-10B, 
and LkHa 336c suggest that the discs of these stars are seen 




h h w 6*110- a*w* a'as* s'o4- ^ha" ter j'ra* 



Figure 9. Continuum subtracted Ho (SHASSA) image of a 
10? 6 X 8? 6 field in the Orion region demonstrating the posi- 
tion of L 1622 (solid ellipse) relative to the association subgroups 
Ori OBla and Ori OBlb (dotted ellipses). Crosses indicate the 
Belt's brightest stars. 



gests that the source of the high pressure which initiated star 
formation in L 1622 may be somewhat more distant than the 
cloud itself. 

The bright, ionized head of the cloud points toward 
Orion's Belt (Fig.O, suggesting that stellar winds and ion- 
ization front from the supergiant members of the associa- 
tion s ubgroup Ori OBlb, whos e age is about 6-7 million 
years (|Blaauwl [l99ll ; iBallvl |2008| ) . might have provided the 
shock which initiated star formation in L 1622. The separa- 
tion of ~ 4.8 deg between L 1622 and the easternmost Belt 
star C Ori corresponds to some 34 pc at a distance of 400 pc. 



nearly edge-on. L 1622-3 exhibits an extremely rich emission 
spectrum. This star might have been born in a small molec- 
ular clump outside the main cloud. IRAS 05522+0146 is a 
Class I source embedded in the dense core L 1622 B. 
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